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We prove that a temperature independent mass distribution is identically zero below a mass 
threshold (mass gap) value, if the pressure satisfies certain inequalities. This supports the finding 
of a minimal mass in quark matter equation of state by numerical estimates and by substitution 
of analytic formulas. We present a few inequalities for the mass distribution based on the Markov 
inequality. 



INTRODUCTION 

Quark matter, searched for in relativistic heavy ion 
collisions, reveals itself in signatures on observed hadron 
spectra which are (can best be) interpreted in terms of 
quark level properties. In particular scaling of the ellip- 
tic flow component vi with the constituent quark content 
of the finally observed mesons and baryons [l[ and suc- 
cessful description of p-r-spectra of pions and antiprotons 
using quark coalescence rules for hadron building [2| uti- 
lize the fast hadronization concept of quark redistribu- 
tion. Albeit this simple idea brings also problems with 
it, e.g. in dealing with energy conservation and entropy 
non-decreasing, these issues can be resolved by using a 
distributed mass quasiparticle model for quark matter 
Q, and are in accord with the quark matter equation 
of state obtained in lattice QCD calculations [4|. The 
surmised mass distribution gives rise to specific equation 
of state (pressure as a function of temperature, p(T)), 
and reversed, a mass distribution may be outlined from 
knowledge on the p(T) curve. 

While traditional, fixed mass quasiparticle models al- 
ready succeed to describe the equation of state obtained 
in lattice QCD [B[ , those mass values are themselves tem- 
perature dependent. Furthermore a temperature depen- 
dent width is associated to the quasiparticle mass, too 
Q. This is unavoidable, as it can be easily seen from the 
following comparison of pressures obtained from a mass 
distribution and from perturbative QCD in the high tem- 
perature limit: 
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mass m defines the scaling function 
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We note that <i>(x) can also be obtained for Bose or 
Fermi distributions instead of the Boltzmann one; the 
numerical difference is overall minor, less than six per 
cent. The basic assumption, (m 2 ) ~ g 2 T 2 sets the 
scale for a simplified treatment of the quark matter pres- 
sure at high temperature. The comparison reveals that 
( m 4 ) 7^ \ rn 2 ) , whence the necessity of a width in 
the mass distribution emerges. The temperature depen- 
dence of the pressure ratio to the massless ideal gas 
value is concentrated on the temperature dependence 
of the coupling constant: g = g(T) in the traditional 
interpretation. We have recently pursued an alterna- 
tive approach to the quasiparticle mass distribution in 
quark matter 0, where a temperature independent 
w(m) distribution is reconstructed from the pressure ra- 
tio cr(T) =p(T)/p id (0,T) curve: 



p CO 

a(T) = J w(m) $ (-) dr 



(3) 



where the ideal gas pressure for relativistic particles with 



This is possible only with a single w(m) curve for each 
exactly known <r(T). Such w(m) distributions show di- 
verging expectation values for positive powers of the mass 
(like ( in 2 )) signaling a high mass tail not falling faster 
than ~ m -3 . Another remarkable property of this ap- 
proach is that it indicates a temperature independent 
threshold (smallest mass) in the w(m) spectrum for lat- 
tice QCD pressure data0,[l|. 

The pressure is, however, not known analytically, the 
numerical results are smeared with error bars. This prob- 
lem is more severe in the light of the fact that eq.® 
is an integral transformation (related to the Meijer K- 
transformation, a generalization of the Laplace transfor- 
mation). There is no mathematical guarantee that by 
inverting such a transformation one obtains close results 
for w(m) from close functions for cr(T). In fact this is 
known as the "inverse imaging problem" [Io| . 
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Our goal in the present paper is to support knowledge 
about a T-independent w(m) mass distribution when the 
pressure p(T) satisfies certain inequalities. In particular 
we prove that if the pressure p(T) is below the corre- 
sponding ideal gas pressure with a given mass Mq at 
all temperatures, then the mass distribution is exactly 
zero for all masses below Mo. For inequalities with other 
than ideal gas pressure curves as estimators we apply 
the Markov inequality for probability measures, which 
directly offers upper bounds on the integrated proba- 
bility P(M) = J Q w(m)dm. It turns out that the ap- 
pearance and value of the highest possible M for which 
P{M) = 0, the mass gap value M , is connected to the 
low temperature behavior of <j(T). Two particular es- 
timators for cr(T), namely $(M /T) with M = 3.2T C 
and exp(— T c /T) are compared to 2+1 flavor lattice QCD 
scaled pressure data in Figure Q] (top) and to pure SU(3) 
lattice gauge theory data (bottom) with M = 2.7T C and 
A = 0.55T C . Of course the temperature scales are differ- 
ent, T c w 165 MeV in the first, T c w 260 MeV in the 
second case. These examples are important for gaining 
a physical insight into the Markov inequality discussed 
below. 
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FIG. 1: Pressure curve estimators for the data from lattice 
QCD simulation of Ref.Q (2+1 flavor QCD) and of [9] (pure 
SU(3) gauge theory). 



MARKOV INEQUALITY AND MASS GAP 

In this section we derive the mass gap based on the 
Markov inequality and use this approach to estimate up- 
per bounds on the integrated probability for the mass 
being lower than a given value. The general form of the 
Markov inequality is given by [ll| 

p({xeX:f(x)>t})<-±r J g(f(x))Mx) (4) 

x£X 

with measure [i, a real valued //-measurable function /, 
and a monotonic growing non-negative measurable real 
function g. The proof, based on the monotonity of in- 
tegration, can be presented in a few lines. For a non- 
negative and monotonic growing function g(t) < g(f(x)) 
for t < f(x). We obtain 



g(t) j d M (x) - / g(t) d/i(a;) <J g ( f(x) ) d»(x). (5) 

f(x)>t f(x)>t f(x)>t 

This quantity can be bounded by 

g(f(x))d»(x) < J g(f(x))d»(x). (6) 

f(x)>t xex 

A division by g(t) > delivers the original statement in 
eq.(0J). 

In order to apply this inequality to the mass spectrum 
we choose f(m) = tM/m. In this case 



and the Markov inequality reads as 
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For a continuum mass spectrum d/i(m) = w(m) dm can 
be chosen with w(m) being the probability density func- 
tion. The generalized Markov inequality stated above 
is valid for general probability measures 12j \x possibly 
including bound state contributions. 



Power law estimate 

In the following we discuss a few examples for mono- 
tonic rising functions g(z), which allow us to draw some 
conclusions about the integrated probability for masses 
below M. Applying the special form of g(t) — t n we 
arrive at 



1 f°° ftM\ n 
P{M) < — - w(m) dm, 



(9) 
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whence we obtain: 



P{M) < M n / m - n w(m) dm. 



(10) 



It is easy to see that the negative integral moments of 
the mass on the right hand side of the above inequality 
are connected to the negative integral moments of scaled 
pressure cr(T) — p(T)/pn(0,T). The final inequality for 
the probability of having masses smaller than M is given 
by 



P(M) < M 
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FIG. 2: Upper bounds for the integrated probability P(M) 
of masses lower than M, based on 2+1 flavor lattice QCD eos 
data0] (top) and for pure SU(3) gauge theory data (bottom). 
Each data line belongs to a fixed value of the scaling param- 
eter t occurring in eq. (|15fl . see legend. In the SU(3) case a 
constant error of 0.02 was assumed in the original p/T 4 data. 



Upper bound by a fixed mass ideal gas 

Let us apply this result to the simplest majorant, that 
of a fixed mass relativistic ideal gas. In this case cr(T) < 
4>(M /T) with some M (cf. dashed line in Figure [TJ. 
Equation (fTTj) leads to 



P(M) < M 



f °° x^^x) dx 



< 



M 
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(12) 



in this case. Should it hold for arbitrary high n, the 
right hand side of this inequality is zero for all M < 
Mq and divergent for M > Mq. In the second case it 
is not restrictive, since P{M) < 1 anyway, in the first 
case this means a mass gap up to Mq. We note that 
this conclusion holds for a general non-negative $(x), for 
which the integrals in eq. (fT2"|) are finite for all n > 0. 
Thus the Bosc-Einstcin or the Fermi-Dirac distribution 
could as well be applied instead of the Boltzmann one. 



Exponential upper bound 

Another possible majorant is the exponential function, 
a(T) < exp(-A/T) (cf. the dotted line in Figure Q] for 
A = T c ). In this case eq. fTTj) delivers 
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The large n limit of this result is given by 
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to leading order in l/n. Again the right hand side ap- 
proaches zero for M < A and diverges for M > A. This 
points out a mass gap stretching to (and including) A 
from zero. 



The pressure itself is a majorant 

The most striking inequality is obtained by using 
g(t) = $(l/t). This function is also admissible, its rise 
from zero to one is strict monotonic. Eq. (fTTj) leads to 



P{M) < 



a(tM) 
WW)' 



(15) 



For t — 1 using the numerical value ^(l) ~ 0.81 one 
arrives at P(M) < 1.23er(M), which can be directly read 
off from numerical simulation or theoretical predictions of 
ct(T). Figure [2] presents curves for different t- values (see 
legend), all being an upper estimate for the integrated 
probability P(M) in the respective cases of 2+1 flavor 
QCD and pure SU(3) gauge theory. The higher seems to 
be the starting Mq value for the rise of the upper bound 
on P(M), the higher also the magnification of the error 
bars. A secure estimate for the P(M) < 0.05 is given for 
masses M > l.7T c = 280 MeV for the 2+1 flavor QCD 
case, while for M > 7.2T C = 1.9 GeV for the pure SU(3) 
gauge case. While in the first case this can be at best 
an avarage between quark and gluon-likc quasiparticlc 
masses, in the second case should be close to observed 
glueball mass. We note that using cr(tM) < $(M /tM) 
in the t — > limit again a mass gap at Mq follows from 
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ea. (fT5|) . In this respect the use of different ojz) functions 
in the Markov inequality does not matter [13j. 

A related version of the inequality (JTSJ) is obtained for 
tM = T, g{z) = §(tM/zT). The upper bound is ob- 
tained at any fixed T as being 



P(M) < 



$(M/T) 



(16) 



Figure[3]plots upper bounds for P(M) obtained using the 
eq. (|16|) . The most restrictive are the lowest temperature 
data for <r(T), they are, however, also the most contam- 
inated by errors. It is probably safe to conclude that as 
much as 90 — 95% of the masses are above 1.5T C w 440 
MeV according to these data. 

Our mathematical treatment of the mass gap leaves the 
point m = in the possible mass distribution as a special 
case. Assuming that there were such a contribution of 
finite measure, i.e. P(0) — a were a finite value between 
zero and one, one concludes from the definition eq.Q 
that in this case a(T) > a would be. There is no sign of 
such an indication in lattice QCD data. 

Finally we note that there is a potential to use our 
method presented in this letter in a context wider than 
quark matter: the quasiparticle test based on the gen- 
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FIG. 3: Upper bounds for the integrated probability of 
masses lower than M based on eq. (|16[) and on 2+1 flavor lat- 
tice QCD eos data0](top) and from Ref. [LPdata] (bottom) 
respectively. The estimates belong to different temperatures 
T near T c indicated in the legend. In the second case a con- 
stant error of 0.01 was assumed in the orig inal p/T 4 data. 



eralized Markov inequality can in principle be done for 
any system with sufficiently known thermal equation of 
state. The estimate for a lowest mass can then be checked 
against knowledge on the mass spectrum obtained from 
the study of correlation functions. 
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a measure normalized to one 

From practical viewpoint, however, in the t — * limit 
the error bars on the origional p/T 4 data are infinitely 
enlarged. 



